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ABSTRACT 



The Herbig Ae/Be star LkHa 101 has been imaged at high angular resolution at a 
number of wavelengths in the near-infrared (from 1 ~ 3 /im) using the Keck 1 Telescope, 
and also observed in the mid-infrared (11.15 //m) using the U.C. Berkeley Infrared Spa- 
tial Interferometer (ISI) . The resolved circular disk with a central hole or cavity reported 
in Tuthill et al. (2001) is confirmed. This is consistent with an almost face-on view (in- 
clination ^35°) onto a luminous pre- or early- main sequence object surrounded by a 
' massive circumstellar disk. With a multiple-epoch study spanning almost four years, 

relative motion of the binary companion has been detected, together with evidence for 
changes in the brightness distribution of the central disk/star. The resolution of the 
Q-f LkHa 101 disk by ISI mid-infrared interferometry constitutes the first such measure- 

ment of a young stellar object in this wavelength region. The angular size was found to 
increase only slowly from 1.6 to 11.15 /im, inconsistent with standard power-law tem- 
perature profiles usually encountered in the literature, supporting instead models with 
a hot inner cavity and relatively rapid transition to a cool or tenuous outer disk. The 
radius of the dust-free inner cavity is consistent with a model of sublimation of dust 
in equilibrium with the stellar radiation field. Measurements from interferometry have 
been combined with published photometry enabling an investigation of the energetics 
and fundamental properties of this prototypical system. 

Subject headings: circumstellar matter - stars: mass loss - techniques: interferometric - stars- 
individual: LkHa 101 
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the mass of the central star (Shu et al. 1990; 
Hollenbach et al. 1994), is sculpted and even- 
tually dissipated by the radiation and wind 
from the newborn star. The structure of these 
disks has been speculative, with uncertainty 
surrounding characteristic sizes and the exis- 
tence of an inner cavity. 

Images of young circumstellar disks in 
Orion(0'dell, Wen & Hu 1993; McCaughrean 
& O'dell 1996) and remnant or fossil disks at 
large distances from the star around "Vega- 
type" sources such as [3 Pic (Smith & Terrile 
1984) have been seen in reflected light, or sil- 
houetted against bright nebulosity. This work 
has attracted intense interest for, as first hy- 
pothesized by Kant (1755), it is likely that our 
own solar system grew out of such a flattened 
primordial nebula (Beckwith & Sargent 1996) 
labelled by Kant an "Urnebel" . 

The framework for understanding stellar 
formation is based on a rapidly-evolving the- 
oretical and observational picture. Although 
the formation of circumstellar disks has long 
been favored, for the case of massive stars 
there have been historical problems with get- 
ting models which fit the shape of the spectral 
energy distribution (SED) (Hillenbrand et al. 
1992) to account for the luminosities in the 
near-IR (Hartmann, Kenyon & Calvet 1993), 
and with forbidden line profiles not matching 
expectations (Bohm &: Catala 1994). It was 
also found that the SED could be fitted by 
spherically symmetric shells (Miroshnichenko, 
Ivezic &; Elitzur 1997; Pezzuto, Strafella & 
Lorenzetti 1997) or composite shell-disk mod- 
els (Miroshnichenko et al. 1999). 

Recent observational evidence has argued 
both for and against the presence of disks. 
Millimeter-wave interferometry has found ev- 
idence for large-scale rotation in Herbig Ae 
stars (Mannings & Sargent 1997, 2000), al- 
beit at much larger spatial scales than the ex- 
pected sizes of the inner disks. However, sim- 
ilar studies of Herbig Be stars (Fuente et al. 



2001) conclude that the disks have dissipated 
at an early stage before the star becomes vis- 
ible. Even restricting attention to just Her- 
big Ae stars, HST coronographic observations 
of AB Aur (Grady et al. 1999) show a more 
circularly symmetric structure consistent with 
disk viewed pole-on, contrasting with the mm 
observations. 

Recent high-resolution imaging in the in- 
frared has also been equivocal. Classical ac- 
cretion disks have generally performed poorly 
in attempting to fit observations from the lat- 
est generation of long-baseline interferometers 
(Millan-Gabet et al. 1999, 2001; Akeson et al. 
2000). More directly, the expected source 
asymmetries from a population of disks, most 
of which must be inclined to the line of sight, 
were not indicated (Millan-Gabet et al. 2001). 
However, of the three systems for which the 
hot inner regions have been well resolved 
by full interferometric imaging in the in- 
frared, HK Tauri B (Koresko 1998), LkHa 101 
(Tuthill et al. 2001) and MWC 349 (Danchi 
et al. 2001), clear disk morphologies (edge-on, 
face-on and edge-on respectively) have been 
established. In contrast to the near-IR where 
disks have proven larger than expected, Hinz 
et al. (2001) failed to detect extended flux 
around a sample of Herbig Ae stars from mid- 
IR nulling interferometry. 

Modellers are now re-thinking the simple 
thin/flared-disk geometries with power-law 
thermal profiles, now seen by many as inade- 
quate in light of the measurements of Millan- 
Gabet and others, and promising new candi- 
dates for reconciling some of the discrepan- 
cies are emerging. One new model (Dulle- 
mond, Dominik, & Natta 2001) involves an 
inner circumstellar cavity whose radius is set 
by thermal dust evaporation which adjoins 
a flaring disk geometry (Chiang & Goldreich 
1997) at larger radii. Such flared disks may 
have thermal instabilities causing runaway lo- 
cal heating, puffed walls and self-shadowing of 
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more distant regions (Dullemond 2000). The 
resultant geometry is capable of intercepting 
more of the stellar radiation, thus boosting 
its near-IR excess while maintaining excellent 
fits to the SED. This has been verified ex- 
perimentally with new spectral data on four 
young stars (Natta et al. 2001). 

LkHa lOlwas first identified as the source 
of illumination of the irregular reflection neb- 
ula NGC 1579 by Herbig (1956), who matched 
the Ha spectrum with a faint, deeply embed- 
ded star lying ~5 arcseconds within the bor- 
der of a dark lane crossing the few-arcminute 
sized nebula. In the years since this discovery, 
it has become one of the most studied young 
stellar objects with observations spanning the 
spectrum. Despite this, pinning down many 
of its basic properties has proved elusive. The 
800 pc distance established by Herbig (1971) 
was based on measurement of two nearby stars 
believed to be in association. Recently, Stine 
& O'Neal (1998) made a strong case based 
on radio photometry that the system is much 
closer; suggesting instead that it is located in 
an extension of the Taurus- Auriga star forma- 
tion complex at 160 pc. 

Part of the difficulty has been the complex- 
ity of the highly anisotropic circumstellar en- 
vironment, with at least four molecular clouds 
(Redman et al. 1986; Barsony et al. 1990) sur- 
rounding LkHa 101 and its associated H II re- 
gion, S222 (Herbig 1956). A deeply-embedded 
population of up to hundreds of probable 
protostars, pre-main-sequence low-mass stars, 
and brown dwarf candidates has been de- 
tected in the radio (Becker & White 1988; 
Stine & O'Neal 1998) and infrared (Barsony, 
Schombert & Kis-Halas 1991; Aspin k, Bar- 
sony 1994). Estimates of the visible extinc- 
tion vary wildly in the literature by almost 
10 magnitudes between values A v = 9.4 (Bar- 
sony et al. 1990) to 18.5 mag (Hou, Jiang & 
Fu 1997). 

This paper presents results from extremely 



high angular resolution imaging obtained with 
interferometric techniques. This has allowed 
observation of the hot, young, self-luminous 
disk surrounding LkHa 101 in the near- and 
mid-infrared. The observational methods and 
the two different instruments used to obtain 
the data are described in Section 2, with the 
results given in Section 3. The astrophysi- 
cal interpretations are discussed in Section 4, 
while section Section 5 contains a summary of 
the important findings. 

2. Observations 

Observations were made at a number of 
epochs over the interval 1997-2001 and uti- 
lized two separate instruments to secure high 
resolution imaging data in the near- and mid- 
infrared. Near-IR observations were taken 
with the Keck-I telescope, while the mid- 
IR interferometry was performed using the 
U.C. Berkeley Infrared Spatial Interferometer 
(ISI). These are briefly described in turn be- 
low. 

2.1. Near-Infrared Interferometry 

Observations at the Keck I telescope em- 
ployed the technique of aperture masking in- 
terferometry in order to recover information 
out to the diffraction limit (Tuthill et al. 2000) 
of the 10 m primary mirror. Placing a mask 
over the telescope pupil has been shown un- 
der some circumstances to confer signal-to- 
noise advantages over full-pupil speckle inter- 
ferometry (Haniff & Buscher 1992). A large 
number of rapid-exposure, high magnification 
data frames were processed to extract Fourier 
amplitudes and closure phases, enabling im- 
ages to be produced from a self-calibration 
algorithm based on the maximum entropy 
method (Sivia 1987; Skilling & Bryan 1984). 
A detailed description of the apparatus and 
data processing can be found in Tuthill et 
al. (2000), and its application in the study of 
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young stars demonstrated (Bittar et al. 2001; 
Tuthill et al. 2001) 

An observing log of near-IR observations 
is given in Table 1 showing the dates, filters, 
and point-source calibrator stars used for each 
data set. The observing wavelength was se- 
lected from the standard complement of inter- 
ference filters in the facility Near-IR Camera 
(NIRC), which offered a range of bandwidths 
(from about 1 to 20%) over the near-infrared 
region (Matthews et al. 1996). The character- 
istics of filters used are given in Table 2. All 
observations utilized the annular ring shaped 
pupil (Tuthill et al. 2000), with the exception 
of 1999 February 05 observation within the 
J-band, which utilized the full unobstructed 
pupil due to the relatively low apparent lu- 
minosity of LkHa 101 in this spectral region. 
An additional difficulty with the J- and some 
H-band observations was associated with the 
fact that the finest fringes were beyond the 
Nyquist sampling limit of the camera. This 
prevented recovery of data at a few of the 
very longest H-band baselines, while for J- 
band even intermediate baselines had serious 
problems with fringe power aliasing, reducing 
the effective resolution to that of a 4m tele- 
scope. 

2.2. Mid-Infrared Interferometry 

Mid-IR visibility data were obtained at 
11.15 fim with the ISI, a two-element het- 
erodyne stellar interferometer located on 
Mt. Wilson, CA. Both telescopes are mounted 
within movable semi-trailers which allowed re- 
configuration of the baseline over the course of 
these measurements. Periodic observations of 
K giant stars a Tau and a Boo were utilized 
to calibrate fringe visibilities and to moni- 
tor drifts to within a few percent accuracy. 
Detailed descriptions of the apparatus, ob- 
serving, and data reduction procedures can 
be found in Hale et al. (2000). A journal of 
LkHa 101 observations taken with the ISI is 



Table 1: Journal of near-IR Keck observations 



Date (UT) 


Filter 


Calibrator Star 


1997 Dec 17 


H 


HD 27349 




KCONT 


HD 27349 




PAHCS 


HD 27349 


1997 Dec 19 


H 


HD 27349 




KCONT 


HD 27349 




PAHCS 


HD 27349 


1998 Sep 30 


H 


54 Per 




BR7 


HD 27349 




CH4 


HD 27349 




PAHCS 


HD 27349 


1999 Jan 06 


H 


54 Per 




CH4 


54 Per 


1999 Feb 05 


J 


HD 28423 


2000 Jan 26 


H 


54 Per 




CH4 


54 Per 


2001 Jul 30 


H 


HD 27349 




CH4 


HD 27349 


2001 Sep 04 


J 


HD 28423 




H 


54 Per 




CH4 


HD 27349 




K 


HD 27349 



provided in Table 3. Spatial frequencies and 
position angles listed are averages of a small 
range spanned during the course of observa- 
tions. 

3. Results 

3.1. Multi- Wavelength Near-Infrared 
Images 

Multi-wavelength images of LkHa 101 are 
given in Figure 1. For the J- and K-band im- 
ages, a single map produced from data taken 
in 2001 Sep is shown. The remaining images 
are averages of four separate observations over 
the interval 1997 December to 1999 January, 
with a stronger weighting given to high signal- 
to-noise data obtained when observing condi- 
tions were favorable. Although this averag- 
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Table 3 

Journal of Mid-IR ISI Observations 



Date 


Tel Separation 


Spatial Frequency 


Position Angle 3, 


Visibility 




(meters) 


(10 5 racT 1 ) 


(Degrees) 




1998 Oct 


4.0 


3.28 


73 


0.964 ± 0.087 


1999 Sep 


13.3 


11.47 


106 


0.648 ± 0.068 



Position angle is measured in degrees East of North 



Table 2: Interference niters used in the near- 
IR 



Name 


Ao 


SX 




w 


(fim) 


J 


1.251 


0.292 


H 


1.658 


0.333 


BR 7 


2.165 


0.022 


K 


2.214 


0.427 


KCONT 


2.260 


0.053 


CH4 


2.269 


0.155 


PAHCS 


3.083 


0.101 



ing of maps could help reduce random noise, 
care had to be taken not to include data from 
widely separate epochs as secular changes in 
appearance have been detected in this source, 
as described in more detail below. 

The main features presented by this system 
are a bright resolved circular disk and a sepa- 
rate unresolved companion about 180 mas to 
the E-NE. The resolved component exhibits 
a central hole or depression in flux, and a 
marked asymmetric brightening towards the 
W-SW limb. These basic features were reli- 
ably reproduced over multiple epochs under a 
range of different experimental conditions. 

It is important, however, to point out in- 
herent imperfections in this imaging process. 
Spurious structures are present in the images, 



for example the faint circular ring appearing 
at the lowest contour level in the K, KCONT 
and CH4 maps at a radius of around 90 mas. 
This feature would be familiar to radio as- 
tronomers 'grating ring' and is caused 
by incomplete Fourier sampling. It should 
also be emphasized that the images are re- 
constructions based on the maximum-entropy 
technique of a target which is not highly 
resolved. There is a very large parameter 
space of images which would fit to the data; 
the maximum-entropy solution attempts to 
choose (in some sense) the smoothest, most 
featureless map from this set (see e.g. Skilling 
& Bryan 1984, Narayan & Nityananda 1986). 
This implicit assumption is not always the 
best one to make in the context of astronomi- 
cal imaging, however for the present we adopt 
it as the standard starting point. The impli- 
cations of using alternate initial assumptions 
are explored in Section 4.4. 

When comparing maps made in different 
colors, it is important to be mindful of changes 
in the performance of the imaging system. 
Moving to longer wavelengths results in the 
loss of system resolution. This is clearly 
demonstrated in the 3/xm (PAHCS) image, 
which is consistent with having the same basic 
morphology seen in the 2 /xm images, however 
the lower resolution is not capable of revealing 
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Fig. 1. — Images of LkHa 101 reconstructed from aperture masking interferometry data taken over 
6 different filter bandpasses. To enhance the signal-to-noise, most frames depicted are the average 
of a number of individual reconstructions taken at various epochs, as described in the text. 
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the fine detail. 

The images at shorter wavelengths, in the J 
and H filters should, on the other hand, have 
higher resolution. Unfortunately, this theoret- 
ical improvement is not realized in the maps 
we see. The main reasons for this are (1) the 
negative impact of the seeing becomes greater 
for shorter wavelengths; (2) the pixel scale of 
the NIRC camera is not adequate to prop- 
erly sample the finest interference fringes in 
this wavelength region; and (3) the target is 
relatively dim, forcing the use of wide band- 
width filters making calibration more difficult 
(Tuthill et al. 2000). In the case of the J filter, 
no mask was used at all resulting in further 
noise from atmospheric turbulence. These 
factors combine to result in a higher level of 
noise in the reconstructed images, and a con- 
sequent loss of detail as compared with the 
K-band maps. 

The set of 2 ^m images (K, KCONT and 
CH4) in Figure 1 present the opportunity to 
examine whether the bandwidth of the obser- 
vation is a factor in determining the appear- 
ance of the image. These filters have simi- 
lar center wavelengths (see Table 2) however 
the fractional bandpasses are approximately 
20%, 7% and 2% for K, KCONT and CH4 
respectively. Some apparent differences were 
found to be due to changes with observing 
epoch (see Section 3.3 below) as the K fil- 
ter image of Figure 1 differs by a couple of 
years from the others. A more detailed quan- 
titative examination, comparing only images 
from the same epochs, revealed no systematic 
differences with filter bandwidth. Indeed, this 
should not be surprising as there are no strong 
spectral features that would imply a different 
origin for light isolated by these three filters. 

3.2. Angular Sizes and Simple Bright- 
ness Models 

Although the images given in Figure 1 do 
give a good representation of the structure 



of LkHa 101, it is desirable to take a step 
backwards towards the raw observables in ex- 
tracting quantitative information. The prime 
reason for this is that fitting model bright- 
ness distributions to the observed visibilities 
and closure phases is a far more direct process 
than measuring intensities and sizes from an 
image reconstructed by a complex nonlinear 
algorithm. 

Furthermore, due to limited Fourier cover- 
age, the mid-infrared interferometry from the 
ISI could not be inverted to produce a map 
for comparison with those in Figure 1. Infor- 
mation on the substantially cooler material to 
which the ISI is sensitive needs to be incorpo- 
rated through a process of the fitting of simple 
model brightness profiles to the data. 

The first set of fits made were to extract 
the position and relative flux of the ~180 mas 
binary companion in the near-IR data. The 
fitting algorithm employed both gradient de- 
scent and grid searching to find the optimum 
chi-squared binary parameters from the cali- 
brated visibility and closure phase data. Rel- 
ative fluxes of the companion are given in the 
second column of Table 4 for wavelengths in 
the K-band and shorter (the companion was 
not detected at long wavelengths) while sepa- 
rations are discussed later in Section 3.3. The 
companion can be seen to be very much bluer 
than the resolved disk component with its rel- 
ative flux rising rapidly from K through to J- 
band. 

Having established the parameters of the 
binary companion, it was possible to artifi- 
cially remove it from the Keck data. This 
was done using the simplistic approximation 
of a binary with two point components, how- 
ever this proved easily adequate for the pur- 
pose of isolating the visibility function of the 
bright resolved disk. Two-dimensional visi- 
bility data were then averaged azimuthally to 
enhance the signal-to-noise. 

Data at four wavelengths through the near- 
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IR (Keck data) and at 11.15 /jm (ISI data) 
are shown in Figure 2. Measurements from 
the entire course of the observing campaign 
have been used, with a higher weighting given 
to the best epochs. Visibilities were fit us- 
ing simple circular Gaussian model intensity 
profiles, with the full- width at half-maximum 
(FWHM) of best fit given in the third column 
of Table 4. Although the images of Figure 1 
show that the intensity profile of LkHa 101 
is better represented by a disk with a central 
depression than a Gaussian, we prefer here a 
simple one-parameter model to give a robust 
estimator of the overall size. Despite its sim- 
plicity, the fits of this simple model in Figure 2 
are quite good. 

The most important thing to notice from 
the Gaussian disk fits of Figure 2 is the dra- 
matically smaller size of the central disk in 
J-band. According to the FWHM sizes given 
in Table 4, LkHa 101 shows a steady slow 
decrease in size over the large wavelength in- 
terval from 11.15 to 1.658 [im, then a sudden 
halving between 1.658 and 1.236 /j,m. Indeed, 
this dramatic reduction is also evident in Fig- 
ure 1 which shows the J-filter image having a 
barely-resolved core compared with images at 
other wavelengths. 

Some increase in size with wavelength is 
to be expected from thermal considerations: 
cooler dust further out will give a greater 
fractional contribution at long wavelengths. 
This probably explains the (~20%) increase 
from 1.658 to 3.083 /an and the (~27%) in- 
crease from 3.083 to the ISI data point at 
11.15 /im. Although thermal radiative trans- 
fer modelling is required to fully exploit the 
implications of these data for the structure 
of LkHa 101, the modest nature of these in- 
creases seems to point to a fairly compact 
bounded structure whose size is therefore not 
a strong function of observing wavelength. 

The most logical explanation for the sud- 
den drop in apparent size from 1.658 and 



1.236 /zm is that the central star itself is be- 
ginning to outshine the disk at shorter wave- 
lengths. Visibility and closure phase data 
severely limit the range of allowed locations 
of this stellar component to near the center 
of the resolved disk. Placed anywhere else it 
would generate strong excursions which were 
not observed. 

Further strong evidence for the presence 
of a hot (blue) point source comes from ex- 
amination of the H-filter visibility curve of 
Figure 2. At high spatial frequencies (above 
~ 4 x 10 6 rad _1 ) the measured data departs 
systematically from the Gaussian disk fit in 
a manner consistent with the presence of a 
weak point source. This is illustrated with the 
dashed curve showing the unresolved compo- 
nent holding up the visibilities at high spatial 
frequency. 

Table 4 shows results from fitting two com- 
ponent "Disk + Point Source" models con- 
sisting of a point-component embedded in a 
Gaussian disk. (Technically, these might be 
regarded as three-component models: recall 
that the binary companion has already been 
subtracted from these data.) This model was 
found to reduce the chi-squared misfit to the 
H-band visibility data by an average factor of 
1.9 over the seven epochs recorded, although 
this improvement comes at the expense of an 
extra free model parameter. 

The preferred model for LkHa 101 can 
therefore be summarized as a resolved disk 
(approximated as a Gaussian of around 40 mas 
FWHM) showing only a modest increase in 
apparent size across the entire near-IR, to- 
gether with a distant blue companion (con- 
tributing 14% at J) and a bright embedded 
blue point source (contributing 55% at J). 
Before proceeding with this interpretation, 
however, it is important to note that alter- 
nate scenarios might be indicated, particu- 
larly with the benefit of full radiative transfer 
modelling and including more complex effects 
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1 .0x1 6 2.0x10 6 3.0X10 6 4.0x10 6 5.0x10 6 6.0x10 6 

Spatial Frequency (rod -1 ) 

Fig. 2. — Visibility curves of LkHa 101 from Keck aperture masking and ISI interferometry data. 
Keck data has had successive visibility offsets of 0.2 added (horizontal dotted lines) in order to 
separate the curves on the page. The 1.236, 1.658 and 2.265 fim curves have had excursions due to 
the presence of the binary companion removed (see text). To enhance the signal-to-noise for the 
Keck data, curves depicted are the average of a number (typically 3-5) of separate observations 
taken at various epochs. Overplotted solid lines show a best-fitting circular Gaussian function to 
the data at each wavelength. For the J and H filter observations, a more complicated 2-component 
model is also shown (dashed line) which is discussed in the text. 
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Table 4 

Simple model fits to visibility data for each of five different observing wavelengths (first column). 
The second column gives the fractional flux of the binary companion, while the third column 
gives the best fitting FWHM of a Gaussian disk model of the central resolved component. The 
last 3 columns give fits of a more complex 2-component model to the J and H-filter data only, as 
described in the text. All errors were computed by observing the statistical spread in values found 
by fitting the relevant models to separate, independently recorded data sets spanning a number of 

epochs. 



Wavelength Binary Companion Gaussian Disk Disk + Point Src 

Relative Flux FWHM Disk Flux Disk FWHM Point Flux 

(//m) (mas) (mas) 



1.236 0.140±.021 16.1±0.8 

I. 658 0.040±.013 36.4±0.8 
2.265 0.008±.003 43.7±1.8 
3.083 45.7±1.0 

II. 15 62.6±10.9 



0.31±0.03 36.8±3.7 0.55±0.03 

0.84±0.02 42.1±2.4 0.12±0.02 



such as scattering which may also play a role 
in altering the appearance of the disk towards 
the blue. 

In a recent imaging study of young mas- 
sive stars, Leinert et al. (2001) reported the 
presence of an extended (0.9") halo contribut- 
ing 0.15 of the flux in H-band. Evidence for 
this was not seen in the visibility curves of 
Figure 2, where a sharp drop would be ex- 
pected over the nearest few points from the 
origin. Interferometric techniques (including 
both the results here and those of Leinert et 
al. 2001) can prove unreliable for the very 
shortest baselines (the so-called "seeing spike 
problem"; e.g. Tuthill et al. 2000). However, 
despite this, a halo as strong as 15% would 
give signals which lie well above this noise pro- 
cess and should have been seen here. We hope 
that discrepancies such as this are sorted out 
as the techniques mature. 



3.3. Time Evolution of Images 

With observations spanning almost four 
years, careful comparisons were made in order 
to determine if LkHa 101 exhibited significant 
changes of morphology with time. In particu- 
lar, the location of the binary companion was 
carefully monitored for evidence of orbital mo- 
tion. A plot of the separation (in right ascen- 
sion and declination) is given in Figure 3. The 
error bars are larger than can be obtained for 
simpler binaries consisting of 2 point-source 
components due to ambiguity introduced by 
the need to locate a fiducial fixed point within 
the complex disk structure. For this purpose, 
the center of the circular disk (ignoring the 
asymmetric limb-brightening) was chosen. 

Figure 3 shows the companion to be mov- 
ing towards the West, with almost all the 
change in the RA coordinate. Fitting a 
uniform velocity to the data yields a mo- 
tion of 6.4±1.3 mas/yr at a position angle 
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Fig. 3. — Relative position of LkHa 101's 
companion as a function of Julian date. Up- 
per panel gives Right Ascension, lower panel 
Declination. The best fit uniform velocity 
model (solid line) is overplotted with the one- 
sigma error in the component velocities indi- 
cated (dashed lines). 



of 275±12 degrees with a starting location 
of 182.5±2.4 mas separation and position an- 
gle 69±2 degrees in 1997 December 17 (JD 
2450799). Data are neither extensive nor pre- 
cise enough to justify fitting of a full orbit 
at this stage. The motion can be directly 
visualized from examination of the images 
presented in Figure 1. The J and K panels, 
whose data were taken in 2001 Sep, show the 
companion visibly further to the West than 
the H, KCONT and CH4 images recorded in 
and around 1998. 

Future monitoring of this companion to de- 
termine the orbit should help in pinning down 
the basic properties of the system, as dis- 
cussed further in Section 4.2. It will also be 
interesting to determine if the companion or- 
bit is coplanar with the disk, as this has im- 
plications for the mechanism of formation of 
the binary system. The interested reader is 
referred to the discussion in Koresko (1998), 
who found HK Tauri B to have a disk which 
was not coplanar with the binary companion. 

Systematic changes of structure were also 
recorded for the resolved disk component of 
LkHa 101. This was studied by compari- 
son of the highest fidelity images, which were 
all produced in the K-band spectral region. 
Images were superior as this represented the 
best compromise between higher theoretical 
resolution (shorter wavelengths), more benign 
atmospheric noise properties (longer wave- 
lengths), and low contamination by the cen- 
tral point-source (longer wavelengths). Fig- 
ure 4 gives a series of K-band images from 
1997 December through 2001 September de- 
picting just the central disk component. 

Although the basic structure of an asym- 
metrically limb-brightened circular disk with 
a central hole or depression is preserved in 
this series, there are more subtle changes 
associated with the peak intensity. It can 
be seen that the brightest quadrant, initially 
to the west-southwest (1997/1998), gradually 



11 




-60-40-20 20 40 60 -60-40-20 20 40 60 
Milliarcseconds Milliarcseconds 



Contours (% of Peak): 1 2 5 10 20 30 50 80 

Fig. 4. — K-band images of the resolved disk of LkHa 101 from data taken at six different epochs. 
With the exception of the KCONT filter which was used in Dec97, all observations were made with 
the CH4 filter. 
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lengthens and shifts to the south (1999/2000), 
before finally splitting into two: a north- 
ern and slightly brighter southern components 
(2001). Although differences as large as this 
might be expected from noisy data in any 
single map, this is unlikely to be the case 
here. These changes are highly systematic 
over time, and multiple contemporaneous ob- 
servations at each individual epoch all corrob- 
orate the average map at that epoch to a high 
degree. 

However, an important caveat on the in- 
terpretation of reconstructed structure near 
the diffraction limit needs to be re-emphasized 
at this point. Although we believe there 
has been real evolution of the morphology of 
LkHa 101's resolved disk, the precise details 
of this change are unclear. The Maximum- 
Entropy picture of the north-south migration 
and splitting of what was formerly a bright 
crescent to the west is only one of many pos- 
sible valid models. This question of ambigu- 
ity over the details of reconstructed images 
will be discussed at greater length in Sec- 
tion 4.4. However despite the uncertainties, 
the changes in the raw visibilities and phases 
which cause the evolution evident in Figure 4 
are real. 

Photometry of LkHa 101 could also be ex- 
tracted from the data utilizing the unresolved 
calibration stars as photometric references. 
Although the experiment was clearly not op- 
timized for such measurements, we were able 
to establish that the infrared fluxes were con- 
stant to within an uncertainty of 10% over the 
course of the observing campaign. Absolute 
flux levels were found to be in good agree- 
ment with published photometry discussed in 
Section 4.3. 



4. Discussion 

4.1. LkHa 101: a face-on accretion 
disk system 

The first qualitative impression given by 
the appearance of the images of LkHa 101 is 
that they conform well to our preconceptions 
for a young system surrounded by an accre- 
tion disk. The interpretation of the circular 
outline of the dust as a face-on view onto a 
disk, rather than a spherical structure, is well 
justified. The inner regions have long been 
known to be highly inhomogeneous as testi- 
fied by the bright reflection nebula NGC 1579 
which must have a low line-of-sight opacity to 
the central star. Further verification comes 
from studies of linewidths and excitation con- 
ditions (Simon & Cassar 1984) revealing the 
existence of a number of regions. A polar 
viewing angle for this star has also been fa- 
vored in the literature (Barsony et al. 1990). 
The degree of circularity was measured by 
fitting to the images with the best signal- 
to-noise, and ignoring the brightest (asym- 
metric) structures. No systematic elongations 
were found in any direction, enabling a limit 
of <35° to be placed on any tilt of the polar 
axis from the line of sight (assuming the disk 
itself to be circular and flat). 

The marked central depression or gap evi- 
dent in the images of Figures 1 and 4 points 
to a hole, or at least optically thin region at 
the center of the disk surrounding the central 
star. Such central cavities have long been in- 
voked in order to help with the fitting of the 
broadband spectra, and are also found in more 
recent models (Hollenbach et al. 1994; Yorke 
& Welz 1996). The latter scenarios apply to 
massive stars upon attainment of their main- 
sequence luminosity, when they commence the 
erosion of their accretion disks in <10 6 yr with 
radiation and wind. It is possible that our im- 
ages of LkHa 101 give the first direct view of 
this process at an early stage, with the star 
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caught in this ephemeral phase having etched 
away the inner disk resulting in the dark hole 
at the center, but with substantial material 
still orbiting in a hot close ring. 

Model C 6=30° A=12 / um 
1 .000 1 I I I I I 



0.667 - 




-1 .000 I I I I I I 

-1.000 -0.667 -0.333 0.000 0.333 0.667 1.000 
S [arcsec] 

Fig. 5. — Model taken from Kessel et al. 
(1998) (Figure 7f) showing a simulated 12 /xm 
image of a photoevaporating disk, viewed at 
an inclination of 30°. 

The final aspect of the image morphology 
which proves to be qualitatively instructive 
is the asymmetric crescent-shaped brighten- 
ing to the southwest, best illustrated in Fig- 
ure 4. The simplest interpretation is that the 
plane of the disk is tilted away from our line 
of sight in the southwest, allowing a clearer 
view to the hottest material on the inner wall 
facing the star in this quadrant. Conversely, 
in the northeast side of the disk, being closer 
to us, cooler material further out is obscuring 
the hot inner parts. Such geometrical argu- 
ments are only consistent with a thick disk, 
torus, or strongly flared structure. Thin flat 
disks, often encountered in the literature, will 
not produce such effects. 

The details of the appearance of a pho- 
toevaporating disk, including a "peculiar 



horseshoe-like feature" for models with an 
inclination of 30°, were predicted in simula- 
tions of Kessel et al. (1998). A model 12 fim 
intensity map extracted from their paper is 
given in Figure 5 showing a strikingly similar 
morphology to our images of LkHa 101. 

The figure was included to illustrate the 
way in which a somewhat inclined, centrally 
illuminated thick torus may make a good 
match to our observations. Although such 
qualitative comparisons are highly encourag- 
ing, thick-disk geometries with inner cavities 
are not unique to photoevaporating models 
and establishing a physical model for this sys- 
tem will require more detailed quantitative 
analysis as discussed below. 

4.2. Basic properties 

Despite extensive observational study, there 
is a great deal of uncertainty over the basic 
properties of LkHa 101. A distance of 800 pc, 
established by Herbig (1971) from observa- 
tions of nearby stars thought to be in asso- 
ciation, had become entrenched in the liter- 
ature. However, this was recently challenged 
by Stine & O'Neal (1998) whose radio pho- 
tometry pointed to a much closer distance, 
suggesting that the star may lie in an ex- 
tension of the Taurus-Auriga star formation 
complex at 160 pc. 

LkHa 101 's bolometric luminosity, esti- 
mated at 4.8 x 10 4 L (Barsony et al. 1990), 
was later revised downwards after it was found 
that up to 75% of the flux might be com- 
ing from a deeply-embedded cluster contain- 
ing hundreds of objects (Barsony, Schombert 
k Kis-Halas 1991). If the 160 pc distance 
scale is adopted instead of 800 pc assumed 
by these earlier authors, the resultant lumi- 
nosity of 4.8 x 1O 2 L0 would imply a mid- 
late B star. Although LkHa 101 may not 
be as highly luminous as the O star in ear- 
lier references, a mid-late B star now seems 
to be an under-estimate. Radio observations 
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showing an ionized shell (Harris 1976; Cohen, 
Bieging & Schwartz 1982), when corrected for 
a distance scale closer than 800 pc, find Ly- 
man continuum photon fluxes as would result 
from a significantly hotter early star (Panagia 
1973). As the photoionizing flux is a strong 
function of spectral type, the radio photom- 
etry seems the most reliable indicator of the 
central star in this case. 

Although revisions to the optical flux such 
as an intermediate distance scale and greater 
relative contribution from LkHa 101 to the 
measured bolometric luminosity might bring 
it more into accord with the radio data, there 
are many difficulties such as known global 
anisotropies and uncertain optical depths 
complicating the picture. However, we pro- 
ceed with the assumption (based on matching 
radio photometry with expected Lyman con- 
tinuum flux from given spectral types) that 
LkHa 101 is an early B0~B0.5 zero age main 
sequence (ZAMS) star which implies an ap- 
proximate mass M* ~ 10 - 20 M ; an effective 
temperature T eff ~ 25 000 - 30 000 K; a lu- 
minosity L ~ 10 000 — 25 000 L ; and a radius 
R* ~ 5i? (Panagia 1973). 

The measured properties of the binary (ap- 
parent separation 180 mas; motion 6.4±1.3 
mas/yr) from Section 3.3 allow us to set 
fully independent constraints on the distance. 
These are derived from the idea that the an- 
gular separation and velocity we see are re- 
lated to the true parameters governing the 
binary motion (linear separation, semi-major 
axis, orbital velocity) by the distance com- 
bined with geometrical foreshortening fac- 
tors. Although we know neither these geo- 
metrical factors nor the true mass, it can be 
seen below that estimates of distance can be 
discriminated by making simple probabilistic 
assumptions about these unknown quantities 
(for example, that we are not looking from an 
unlikely, highly foreshortened perspective at 
the binary). 



Firstly if we take the 800 pc distance scale, 
then the orbital period of a binary with semi- 
major axis 180 mas and an assumed system 
mass of 15 M is 446 years. On the plane 
of the sky this would appear as a motion of 
2 . 5 mas /yr . There number of ways our 

measurement of 6.4 mas/yr could easily un- 
derestimate the true motion, however it is 
more difficult to see how we have found a 
larger apparent motion. The enclosed system 
masses would need to be 100 M Q in order to 
generate the observed apparent motion at this 
distance. The companion could be approach- 
ing periastron in a highly eccentric orbit and 
thereby moving faster than expected for its 
separation, or alternatively the objects may 
not form a bound pair at all in which case 
the apparent superposition we see will last 
<100yr. Although they cannot be ruled out, 
none of these explanations seem very satisfac- 
tory and we conclude LkHa 101 is likely closer 
than 800 pc. 

If we consider now the 160 pc distance, an 
apparent motion of 28 mas/yr (orbital period 
40 yr) would be seen for a 15 M Q circular 
binary in the plane of the sky with radius 
180 mas. As mentioned above, there are a 
number of ways in which our measurements 
may underestimate the true state of affairs. 
Apparent motions as low as 6.4 mas/yr could 
be due to a large hidden line-of-sight compo- 
nent to the velocity vector, although this is 
somewhat unlikely as it requires close align- 
ment of the present orbital direction to the 
line-of-sight. The separation of 180 mas used 
for the orbital computation might similarly be 
an underestimate due to foreshortening (the 
companion has a much larger separation, but 
on a vector close to our line-of-sight), however 
this is statistically extremely unlikely, requir- 
ing much closer alignments. Allowing system 
masses as low as 5 and 10 Mq results in pe- 
riods of 69 and 49 yr giving apparent motions 
of 16 and 23 mas/yr. The 5M assumption 
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can match the observations without highly- 
unlikely inclinations, however it is doubtful 
that stars at this light end of this range are 
energetic enough to power the observed HII 
region. As before, it is also possible to fit the 
observations by assuming a highly eccentric 
orbit: in this case the binary would need to 
have been observed moving slowly near apas- 
tron. 

If we instead invert the question and ask 
what is the distance given the apparent mo- 
tion, we arrive at 430 pc for a 15 M & circu- 
lar binary orbiting in the plane of the sky. 
However, we can do a little better than this, 
at least in a statistical sense, by making 
the assumption that we have observed one- 
dimensional projections of velocity and sepa- 
ration vectors taken from a random isotropic 
population (resulting in likely underestima- 
tion by a factor of 7r/4). Accounting for this 
results in a likely distance of ~ 340 pc, al- 
though values in the range 200-500 pc can 
be comfortably accommodated given the un- 
certainties in the mass and projection angle, 
and the experimental error on the apparent 
motion. Allowing for eccentricity widens this 
range still further, although numerical simula- 
tions fitting for various orbits did find far more 
solutions at the intermediate distance than at 
the extremes. Unfortunately, the 180 yr pe- 
riod of such a system makes it unlikely that 
precise orbital elements, lifting some of the 
uncertainty in these calculations, can be de- 
rived soon. 

In conclusion, although we can rule out 
neither of the earlier distance estimates, 
our imaging results favor an intermediate 
~ 340 pc scale. As the previous distance 
scales have not been established with any de- 
gree of certainty, we include this new value 
despite the large uncertainties in its deriva- 
tion above. For the remainder of this paper, 
interpretations based on the short (Stine & 
O'Neal 1998), intermediate (this work) and 



long (Herbig 1971) scales are given. 

4.3. Spectral Energy Distribution 

The observed spectral energy distribution 
for LkHa 101 is given in Figure 6, where the 
data have been taken from various literature 
sources. The overplotted 750 K blackbody 
curve was found (Cohen & Woolf 1971; Da- 
nen, Gwinn, & Bloemhof 1995) to fit the near 
to mid-infrared spectral regions. We were able 
to extract photometry from our own obser- 
vations, and found the fluxes to be generally 
consistent with the literature values, in partic- 
ular good agreement was found with the near- 
infrared data of Bergner et al. (1995). We 
therefore presume that these historical data 
are fairly representative of the present epochs 
under study. 

Wavelength (/Am) 
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Fig. 6. — Spectral energy distribution of 
LkHa 101 from the optical through the mid- 
infrared taken from various literature sources 
(as indicated). The overplotted solid line is a 
750 K blackbody curve (Cohen & Woolf 1971; 
Danen, Gwinn, &: Bloemhof 1995) 

An analysis of the energetics of LkHa 101 
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and its associated region is complicated by the 
fact that it is known to lie behind thick clouds 
of obscuring material. The spectrum must 
therefore be corrected for the wavelength- 
dependent extinction (reddening). This has 
been done using dust constants appropriate 
for so-called "outer cloud dust" (Mathis 1990) 
expected in regions of heavy obscuration. 

Widely disparate estimates of the line-of- 
sight visual extinction (A v ) can be found in 
the literature, with values of 9.4, 9.7, 11.2, 
14.2, 15.7, Sz 18.5 obtained from, respectively, 
Barsony et al. (1990); McGregor, Persson, 
& Cohen (1984); Kelly, Rieke & Campbell 
(1994); Thompson et al. (1977); Rudy, Er- 
win, Rossano, & Puetter (1991); Hou, Jiang 
& Fu (1997). In the following discussion, we 
proceed with the relatively recent, intermedi- 
ate value of 11.2 from Kelly, Rieke & Camp- 
bell (1994), although extinctions encompass- 
ing the entire range were trialed. When used 
to de-redden the spectrum of LkHa 101, val- 
ues to the upper end of the range were found 
to give unrealistic spectra rising very sharply 
in the blue. At the other extreme (low A v ) 
the system could be modelled, but a relatively 
small, cool central star was found making it 
difficult to reconcile with the radio measure- 
ments of Lyman continuum photon fluxes. 

The de-reddened spectrum of LkHa 101, 
showing only the visible and near-infrared, is 
shown in Figure 7. The spectral data points 
(diamonds) were taken from Figure 6, averag- 
ing multiple measurements at the same wave- 
length where necessary. There is now a strong 
indication of a 2-component curve to the SED, 
with a red spectrum dominating wavelengths 
longer than K band, while visible regions ap- 
pear to betray the presence of a hot blue ob- 
ject. The transition region between these two 
cases falls around 1 fim. 

Utilizing information from the imaging and 
model fitting in Section 3, we are able to de- 
compose the full SED of Figure 7 into sep- 



arate component spectra, each corresponding 
to individual features identified using the high 
resolution techniques. Firstly, we are able to 
extract the spectrum of the binary compan- 
ion by noting that we have obtained measure- 
ments of its relative flux in the J, H and K 
bands as given in Table 4. The resultant SED 
of the companion, isolated in this fashion, is 
plotted as triangular points. 

In a similar fashion, the remaining flux 
from the total can be divided into compo- 
nents due to the disk and central point source, 
where the flux ratios have again been taken 
from Table 4. The central point source, orig- 
inally identified in Section 3.2 from Figure 2, 
is indicated by stars, while the resolved com- 
ponent of the flux assumed to come from the 
disk is indicated by open circles. Note that 
the central point component was only clearly 
identified at 2 wavelengths, J and H, but we 
have also overplotted a central point compo- 
nent which contributes 2.5% of the flux at K 
band. Although there may be weak evidence 
for such a signal in the visibility curve of Fig- 
ure 2, this was done mainly as an illustration 
of the expected K-band flux of the central 
component extrapolated from J and H. For 
wavelengths longer than K, the circumstellar 
disk is expected to dominate the SED. 

For each of the three components of the 
system identified, a blackbody curve has been 
plotted. Both stellar components have been 
modelled as having a temperature of 25 000 K, 
(B0. 5 ZAMS star from Panagia 1973) al- 
though the limited spectral coverage and un- 
certainties on the photometry mean that this 
is not well constrained by the data. A near- 
infrared dust temperature of 1000 K has been 
used to fit the resolved disk component, and 
in contrast to the stellar case, the temper- 
ature is fairly well constrained to lie within 
±50 K of this value. The model spectral en- 
ergy distributions for the central star, disk, 
and binary companion are given in Figure 7 
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Fig. 7. — Spectral energy distribution of LkHa 101 corrected for wavelength-dependent extinction 
with A v = 11.2 (Kelly, Rieke & Campbell 1994). Using information from the imaging, the total SED 
(diamonds) has been divided among the components of the LkHa 101 system: the binary companion 
(triangles), disk (open circles) and central point source (stars). These separate components have 
also each been modelled by blackbody spectra, as overplotted in the figure (see text for details). 
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in dotted, dashed and dot-dashed linetypes 
respectively. 

Constraining the angular sizes to match 
the recorded flux levels, we deduce that the 
dust shell must subtend an apparent area of 
800 mas 2 on the sky (an area equivalent to a 
uniformly illuminated circle of radius 16 mas). 
This is in agreement, to within a factor of 
~ 2, with the flux- weighted observed angu- 
lar size of the dust cloud (~ 1500 mas 2 ; ra- 
dius of circle of brightest emission 21 mas) 
from the near-infrared imaging, and to within 
a greater factor (~ 4) from the mid-infrared 
size. These discrepancies could probably be 
taken into account by more complex models 
in which a variation in disk temperature and 
opacity with radius is allowed. 

The central star and binary companion 
have apparent diameters of 56 and 29 fias (as- 
suming the 25000 K temperature). Within the 
context of this model, we are able to analyze 
the energetics of the system, finding about 
~ 20% of the bolometric flux is reprocessed 
radiation from the dusty disk, while the bi- 
nary companion contributes some ~ 15%. It 
is important to emphasize that these conclu- 
sions are dependent upon a number of prior 
assumptions, for example the visible extinc- 
tion and the temperatures of the stars (in 
particular, the temperature of the compan- 
ion is all but unknown). Although we have 
taken reasonable estimates for these quanti- 
ties, the high degree of uncertainty means 
that the resulting model should be regarded 
clS Si plausible exploration rather than an op- 
timized global fit. 

It is now possible to use the prior estimates 
of stellar distance to translate from apparent 
to absolute quantities. This has been done in 
Table 5 for the three different distance scales 
discussed in Section 4.2. The absolute lumi- 
nosities (expressed in 10 3 L Q ) and radii of the 
central source, the binary companion, and cir- 
cumstellar disk are given on consecutive pairs 



of rows. It is very interesting to note that our 
340 pc distance scale results in a stellar lumi- 
nosity and radius appropriate for a Bl ZAMS 
star (Panagia 1973): in accord with our expec- 
tations from radio photometry, and quite close 
to (although slightly smaller than) our ini- 
tial model assumption (B0. 5; L = 11000 L Q ; 
R = 5.1 Rq). Note that neither the shorter 
nor the longer distance scales give this re- 
assuring internal consistency for the model, 
however as noted above, different choices of 
A v or stellar temperature might plausibly re- 
sult in models which work at these distances. 

The value of Rdust given is the diameter 
of a 1000 K circular disk which would gener- 
ate the observed blackbody flux (note this is 
a simple circular area on the sky). R ca vity 
gives the expected linear size of the observed 
21 mas radius circular ridge of brightest emis- 
sion from the images of Figure 4. The close 
agreement between the observed and expected 
sizes (regardless of distance) gives confidence 
in the modeling of the SED and our choice of 
A v . It should be emphasized, however, that 
this is a very simple blackbody model and a 
more detailed treatment should try to find a 
self-consistent thermal profile of the disk. 

We identify the circular structure as the 
hot inner walls of the cavity whose size is 
likely set primarily by direct radiative dust 
sublimation from the central star (e.g. Tuthill 
et al. 2001) rather than disk-reprocessing 
or viscous-heating processes usually assumed 
(Lynden-Bell & Pringle 1974; Hillenbrand et 
al. 1992). Finally, Table 5 gives R su blimation, 
the 1500 K sublimation radius of isolated dust 
grains for the system under the assumptions 
implicit for the 3 distance scales. Grains sur- 
viving at the inner edge would be those best 
able to reradiate the stellar field, so the as- 
sumption of blackbody grains has been made. 

It can be seen from the last two lines in 
the table that the sublimation radius is well 
matched to the size of the inner cavity, irre- 
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Table 5 

Derived physical parameters for the LkHa 101 system, computed for 3 different possible stellar 

distances. 



Distance 


160 pc 


340 pc 


800 pc 


L'star 


1300 L 


5900 L© 


30000 L Q 


Rstar 


1.9 R Q 


4.1 R Q 


9.7 R Q 


L companion 


340 L Q 


1500 L© 


8500 L© 


Rcompanion 


1.0 R e 


2.1 Rq 


5.0 R Q 


Ldust 


270 L© 


1200 L© 


6700 L© 


Rdust 


2.6 AU 


5.4 AU 


12.7 AU 


Rcavity 


3.4 AU 


7.1 AU 


16.8 AU 


Rsublimation 


2.5 AU 


5.3 AU 


12.5 AU 



spective of the assumed distance. This agree- 
ment is encouraging for a simple radiatively 
set inner cavity model, although clearly a 
more detailed treatment could also incorpo- 
rate effects such as disk self-radiation (Bell 
1999). It is interesting to note that mate- 
rial must be continually evaporated at the in- 
ner edge in order to maintain this radius; this 
gas is presumably either ultimately accreted 
by the central star or swept out into the sur- 
rounding HII region. 

4.4. Parametric Imaging 

As we touched upon in earlier sections, the 
images produced by the Maximum Entropy 
method are not unique, nor do the underlying 
algorithms always represent the best possible 
approach. In particular, it may be possible to 
bring extra information to the mapping pro- 
cess, rather than making the assumption that 
the region of sky to be imaged is most likely 
flat and featureless. 

To take a concrete example, we can in- 
fer from Figure 7 that the images in the ri- 
band should have a point source at the center 



contributing some ~ 2.5% of the flux. Even 
concentrated into a single diffraction-limited 
beam, a feature as weak as this could not be 
seen amidst the bright structure within the 
resolved disk. This is verified upon exami- 
nation of images such as those in Figure 4, 
where no trace of the central point-source is 
seen. However, even bright point sources can 
fail to appear in MEM maps as the algorithm 
will penalize fine structure in the image if al- 
ternate smoother maps can be found which 
also fit the data. 

In this case, we are able to bring our knowl- 
edge of the likely presence of a point source 
into the mapping process. This is done by us- 
ing a prior, or default map (e.g. Monnier et 
al. 2001). The prior is the map the algorithm 
will default to in the absence of any constraint 
from data. We have used a prior map which 
contains our expected 2.5% point source at 
the center, thereby permitting MEM to add 
fine structure without penalty. The resultant 
image is given in Figure 8. 

The appearance of a compact central com- 
ponent, in accord with our physical expecta- 
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Fig. 8. — Image made with a 2.5% point- 
source prior from data taken in 1998 Septem- 
ber through the CH4 filter. 

tions for this system, teaches a valuable lesson 
in caution when interpreting maps made from 
data taken at close to the diffraction limit. 
For partially-resolved targets, images cannot 
be definitive but have some model-dependent 
structure. However, this is not to say that 
the images are in any sense arbitrary. A good 
example is the almost precisely central loca- 
tion of the point-source in the Figure 8 im- 
age. Although many trials with priors of var- 
ious sorts were tried, all cases where there 
was a strong point-source in the output map 
had it located at the center. This not only 
agrees with our physical expectations for the 
system, it is also required by the data. An 
off-center point source would generate strong 
telltale visibility and phase signals as was the 
case for the 180 mas binary already mapped. 

Features for which there is strong evidence 
in the data will always appear. An example is 
the limb-brightened ring, which was present, 
almost unchanged, for all images made. More 
subtle things, such as the changes in the lo- 
cation of the brightest part of the disk seen 
in Figure 4, are more open to interpretation. 
It is apparent that some form of systematic 
change in the brightness distribution is un- 



derway. As the orbital period of material in 
the disk is about 5 yr assuming the 340 pc dis- 
tance, it is tempting to interpret the changes 
as being due to dynamical evolution of the 
disk. Features such as spiral density pertur- 
bations were predicted by Adams, Ruden, & 
Shu (1989), and such instabilities may be aug- 
mented by the presence of a binary compan- 
ion, or accounting for self-gravity in a massive 
disk (e.g. Woodward et al. 1994). 

However, the actual depiction in the figure 
of the bright western crescent splitting to the 
north and south is model-dependent. For ex- 
ample, the location of the brightest quadrant 
on the disk was found to be affected by small 
variations in the location and brightness of 
the central star. This leads to the possibility 
that things unrelated to the disk itself, such as 
varying line-of-sight extinction, may generate 
changes in appearance of the reconstructed 
image. We refrain, therefore, from specula- 
tion on the detailed cause of these changes. 

5. Conclusions 

Near-infrared images of the Herbig Ae/Be 
star LkHa 101 have been obtained from a 
multi-epoch study utilizing interferometry on 
the Keck 1 telescope to obtain information at 
the diffraction limit (tens of milli-arcseconds). 
The mid-infrared U.C. Berkeley Infrared Spa- 
tial Interferometer was also able to resolve the 
system at 11.15 //m. LkHa 101 presents a re- 
solved circular disk with a central hole or cav- 
ity. A relatively blue-spectrum binary com- 
panion 180 mas to the E-NE is also reported. 

The morphology of the resolved limb- 
brightened ring is interpreted close to 
face-on viewing angle onto a geometrically 
thick torus or circumstellar disk with a hot 
inner wall facing the central star. The size 
of the cavity is consistent with the radiative 
equilibrium temperature for dust sublimation. 
A relatively slow increase in apparent diam- 



21 



eter of the disk with observing wavelength 
over a decade from 1.2 to 11.15 /im implies 
circumstellar material with a relatively com- 
pact density and/or thermal profile, arguing 
against classical power-law temperature pro- 
files usually encountered in the literature. 

Relative motion of the binary compan- 
ion, together with simple assumptions on 
the masses and geometrical projection, fa- 
vor a likely intermediate distance of around 
~340pc, although the earlier large (800 pc) 
and particularly small (160 pc) suggested dis- 
tance scales cannot be ruled out. 

Combining the interferometry with pub- 
lished photometry enables a decomposition of 
the spectral energy distribution into 3 compo- 
nent parts in the infrared. In addition to the 
180 mas companion and the resolved disk, we 
report the probable isolation of a bright unre- 
solved source at the center of the cavity in the 
disk which has been identified as the primary. 
This is likely an early B star and the source 
of the photoionizing radiation driving the HII 
region. With more detailed radiative transfer 
modelling and studies at yet higher resolution 
from the latest generation of optical interfer- 
ometers, this fascinating system is among the 
most promising candidates for in-depth study 
of the workings of a massive young star. 
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